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Abstract 

We calculated the next-to-leading order (NLO) QCD corrections to the cross sections for direct 
top quark productions induced by model-independent flavor changing neutral current couplings at 
hadron colliders. The NLO results increase the experimental sensitivity to the anomalous couplings. 
Our results show that the NLO QCD corrections enhance the leading order (LO) total cross sections 
at the Tevatron Run 2 about 60% for both of k+ c and k+ u couplings, and enhance the LO total 
cross sections at the LHC about 40% for k\ c couplings and 50% for k\ u couplings, respectively. 
Moreover, the NLO QCD corrections vastly reduce the dependence of the total cross sections on 
the renormalization or factorization scale, which leads to increased confidence in predictions based 
on these results. 
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I. INTRODUCTION 



The flavor dynamics, such as the mixing of three generation fermions and their large 
mass differences observed, is still a great mystery in particle physics today. The heaviest 
one of three generation fermions is the top quark with the mass close to the electroweak 
(EW) symmetry breaking scale. Therefore, it can play the role of a wonderful probe for 
the EW breaking mechanism and new physics beyond the standard model (SM) through 
its decays and productions. An important aspect of the top quark physics is to investigate 
anomalous flavor changing neutral current (FCNC) couplings. Within the SM, the FCNC 
processes are forbidden at the tree-level and highly suppressed by the Glashow-Iliopoulos- 
Maiani mechanism at one-loop level. All of the current precise measurements of various 
FCNC processes, for example, b —> sj |2j, agree with the SM predictions well. However, it 
should be noted that only few experimental constraints on the top quark FCNC processes 
are available so far, and many new physics models allow the existence of the tree-level FCNC 
couplings, and may greatly enhance some FCNC processes. Since the CERN Large Hadron 
Collider (LHC) can produce abundant top quark events, the measurement of the top quark 
rare processes will become possible. Therefore, a good understanding of the theoretical 
predictions of top quark FCNC processes is important for searching new physics. 

Any new physics effect involved in top quark FCNC processes can be incorporated into 
an effective Lagrangian in a model independent way 0, 0] : 

q=u,c q=u,c 

where A is the new physics scale, 9w is the Weinberg angle, and T a are the Gell-Mann 
matrices satisfying Tr (T a T b ) = 5 ab /2. The vf q ,af q can be complex in general, and vf q = 
af q = in the SM. = d^G a u - d u G* - g s ] ahc G\G c v and similarly for the other gauge 
bosons, k is normalized to be real and positive and /, h to be complex numbers satisfying 
for each term |/| 2 + |/i| 2 = l. 

The top quark FCNC processes induced by various anomalous couplings have been studied 
in detail (for a most recent review, see 0). In general, top quark decay processes provide 
the best place to discover top FCNC interactions involving anomalous t — q — 7 and t — q — Z 
couplings. However, for t — q — g anomalous couplings, i.e, the forth term in Eq.(l), the best 



processes are the top quark FCNC productions (6J, |2j , and the direct top quark production 
processes are the most sensitive ones Q. The current constraints on the t — q — g anomalous 
couplings are < 0.47 TeV -1 Q, and as pointed out in |4j, |2| based on tree-level analysis 
of the direct top quark productions, the couplings can further be detected down to — 
0.0084 TeV" 1 , = 0.0033 TeV" 1 at the LHC, and & = 0.062 TeV" 1 , & = 0.019 TeV" 1 at 
the Tevatron Run 2. However, since these processes involve strongly interacting particles in 
their initial states and final states, the lowest order results in general have large theoretical 
uncertainties. In this paper, we present the complete NLO QCD calculation for the cross 
section of direct top quark productions induced by model-independent FCNC couplings at 
hadron colliders. 

The arrangement of this paper is as follows. In Sec. II we show the LO results. In Sec. Ill 
we present the details of the NLO calculations including the virtual and real corrections. In 
Sec. IV by a numerical analysis we present the predictions for total cross sections at hadron 
colliders. 

II. THE LEADING ORDER RESULTS 

We start from the flavor- changing operators including t — q — g anomalous couplings: 

- 9s E ^ta^T\f tq + %h\ ql ,)qG% + H.c, (2) 

q=u,c 

and consider the leading order (LO) scattering process of hadrons P\ and P 2 to top quark: 
Pi(pi) +P2G02) -+g(h) +q(k 2 ) -+t(p t ) +X , q = coiu. (3) 
The spin and color averaged squared matrix element is 

l-^ B | = -3- \ -f J m u (4) 

which, after the phase space integration, gives the LO partonic cross section (the masses of 
charm quark and up quark are neglected) 

o\z) = ^{^ m \8il-z\ (5) 

where s = (ki + k 2 ) 2 and z = mf/s. The LO total cross section at hadron colliders is 
obtained by convoluting the partonic cross section with the parton distribution functions 
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(PDFs) Gg tq /p in the proton(anti-proton): 

o- B = J dx 1 dx 2 [G g / Pl (xi, /jL F )G q/ p 2 (x2, hf) + (xi a; 2 )](T B . (6) 
where fip is the factorization scale. 

III. THE NEXT-TO-LEADING ORDER CALCULATIONS 

For convenience, in this section we present only the results for gc —>■ t, and the ones 
for gu — > t can be obtained by replacing c with u in related formulas. The cross sections 
beyond LO involve the computation of one loop virtual gluon corrections [Figs^Vi)-^t> 12)] 
and real gluon bremsstrahlung contributions to leading order processes [FigsJ^ri)-^^)]. 
We also include processes with two gluons in the initial states [FigsJ^Pi)-^^)] as well as 
two quarks in the initial states (Figs|2f si)H2ts3)). We will encounter ultraviolet (UV) and 
infrared (IR) (soft and collinear) divergences in our computation. We have used n = 4 — 2e 
dimensional regularization [9] to regulate both these divergences, and all divergences appear 
as l/e a with a = 1, 2. 

The virtual corrections to direct top productions consist of self-energy and vertex di- 
agrams, which represent the QCD corrections arising from quarks and gluons. The total 
virtual amplitude contains UV and IR divergences (the imaginary part is neglected): 

= -^oY^-- + — + — -15W 

12vr \ ej R e IR e uv 3 J 

+ (ISZ® + Uzf + Uzf + 8Z gs + 5Z k9j ^ M b , (7) 

where C e = T(l + e) (-^r 1 ) , and the UV divergences are renormalized by introducing coun- 
terterms for the wave function of the external fields (SZ 2 ,8Z 2 , 5Z 2 ), and the coupling 
constants (5Z gg ,SZ K 9j A ). We define these counterterms according to the following conven- 
tion. For the external fields, we fix 5z!f \ SZ^ and SZ^ using on-shell subtraction, and, 
therefore, they also have IR singularities: 

= -^C e ( n -L --)(—-—)- , 
2w \ 3 2/ \euv ^-m' 67r e uv 
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37r \e uv e IR 

SZ? = -^Cj— + — + A). (8) 

37T \e uv e IR 



where n/ = 5. For the renormalization of g s , we use the MS scheme modified to decouple 
the top quark ^(|, i.e. the first rif light flavors are subtracted using the MS scheme, while 
the divergences associated with the top-quark loop are subtracted at zero momentum: 

5Z q = ^ir(l + e )(4vr) e (^-— )— + -^C e — . (9) 
93 4tt V A ; V 3 2 ! e uv 12tt e uv V ; 

Thus, in this scheme, the renormalized strong coupling constant a s evolves with n/ light 
flavors. Finally, for 6Z K gj A , we adopt the MS scheme and adjust it to cancel the remaining 
UV divergences exactly: 

at 1 
£ir(l + e)(47r) e — . 
07T e uv 

Therefore, now A4 vnt has no UV-singular contributions. 

The partonic cross section of real gluon bremsstrahlung are 



sz *j* = ^r(i + £ )(4*)'— . (io) 
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+ — )S(1 - z) + — [I8z 2 - Uz + 40 - - 
^ e IR e IRJ e IR\ z 



—^- \ + 85(1 -z) + -T-. 1 r f 27z 2 - 57z - 2{llz + 27) In z - 23 



11 \ 77z + 27 An(l -z)\ ( , 83 36 \ , 



/no 2 on 107 18 2 V , , a 73^ 63 109 
+ 28z 2 - 22,2 H In z + 14z 2 - 



,(H) 



4 z 1-*/ 8 4 £ 

where C' e = §nE^j (^ffi) ■ The "plus" functions appearing in the above results are the 



1 i; I ' ^ 

r(l-2e) 

distributions which satisfy the following equation: 



dzf + (z)g(z)= f dzf(z)(g(z)-g(l)), (12) 



-1 rl 

, dz f+(z) g(z) = \ 
/o J 

where g(z) is any well-behaved function in the region < z < 1. 

Combining the contributions of the LO result, the virtual corrections and the real gluon 
bremsstrahlung, we obtain the bare NLO partonic cross section: 

-bare/ -i / \ 87T 2 a s ( «f c \ 2 r/-i \ 

Vgc = -^—[^] rn t S{l-z) 



+ 



4vr 2 



- 15 



29 \ ix\ 



6(1 -z) + 



5(1 - 2) 



27z 2 -57z-2(77z + 27) lnz 



-23- 



ll\ 77z + 27/ln(l 



V 1 ~ z 



36^-28* + —- — ) ln(l-z) 
2 z 



+ 28r - 22z + 



107 1? 



2 \ , 2 73z 63 109 

1 In ^ + 14z 2 - ■ 



+ -r ~ 



; (13) 



4 z 1 - zj 8 4 8z y 

where N c = 3. Now the soft divergences coming from virtual gluons and bremsstrahlung 



contributions have cancelled exactly according to the Bloch-Nordsieck theorem [111 ]. The 
remaining divergences are collinear. 

Moreover, calculating the processes with two gluons or two light quarks in the initial 
states, we obtain the following partonic cross sections 

1 



lira 2 2 

— — m: 

3s A * 



, z* + (l-zY 1 (± V[z 

2 e IR + \A + 8 4(1 + *) 



3 \ i /o 2 n iM , 259 2 505 1 5 

+- In z + (2z 2 - 2z + 1) ln(l - z) z 2 H z H 1 

8/ V M ; 64 128 16 128z 



<7^ e (z, 1/e^) 



327mJ/^ c y 2 

— — m t 

9s A * 



1 + (I - z) 2 1 
-C f 1 ' — + 2zln(l-z) -zlnz 



10 4 2 7 
-4 ln(l - z) + 2 In z H h - ln(l - z) In z 

3 z z 3z 



a c b r(z,l/ e/R ) 



167ra 2 / k 9 ^ 



9s 



tc \ 2 

A ^ 



l + Q-z) 2 l 

-C'— ^= ^— + 2zln(l-z) -zlnz 

z e IR 



2 4 7 4 2 2 
— z - 4 ln(l - z) + 2 lnz 1 h - ln(l - z) In z + -z 2 

3 V ; 3z 3 z V ; z 3 



*5&(MA/ii) 



167ra 2 / ~ 
9s V A 



m. 



-C 



+ 2z ln(l — z) — zlnz 



4 2 2 

-41n(l - z) + 21nz + 3 + -ln(l - z) - - lnz - - 

z z z 



9s I A 



,'2 

,2 / z 2 



m f \ —z — z 
' \3 



1 

3z~ 



(14) 
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Note that the divergences appearing in above expressions are collinear, and the patonic cross 
section of the qq(q ^ c) initial state is free of singularities. 

The bare partonic cross sections in Eqs. and (J 14)) . which contain the collinear singu- 
larities generated by the radiation of gluons and massless quarks, have a universal structure, 
and can be factorized into the following form to all orders of perturbation theory: 

(j^ re (z, l/em) = T ca(z, hf, l/e IR ) <g> T db (z, fx F , l/e IR ) <g> a cd (z, /i F ) , (15) 

c,d 

where [i F is the factorization scale and ® is the convolution symbol defined as 

f{z)®g{z) = f ^-f(y)g( z -)- (16) 
J z y \y ) 

The universal splitting functions T cd (z,fip, 1/em) represent the probability of finding a par- 
ticle c with fraction z of the longitudinal momentum inside the parent particle d at the scale 
Hp. They contain the collinear divergences, and can be absorbed into the redefinition of the 



PDF according to mass factorization 
have to 0(a s ) 

T^z^F^/em) = S cd S(l-z) 



12j. Adopting the MS mass-factorization scheme, we 



1 as r(i-e) 

e IR 2tt T(l - 2e) 



(17) 



where P^\z) are the leading order Altarelli-Parisi splitting functions 



l + z2 3 r/, n 



[l-zf + z 2 



4(l-z) 2 + l 
3 z 



Pgg( k 



2N r 



1 - z 



l-z) 



+ z(l-z) 



+ (y-f)^(l-^). (18) 



After absorbing the splitting functions T cd (z, fip, 1/^ir) into the redefinition of the PDFs 
through the mass factorization in this way, we have the hard-scattering cross sections 
a a b(z,fi F ), which are free of collinear divergences, and depend on the scale fip'- 



57T a* / K 



a gc {z,/j,F) 



3s 



A" 



m^8(l — z) 
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2N r 



\l-z) + z 



f - — - + z{l - z) 



+ {T-i) s(1 - z)+ 3{w^r + + 2 S{1 - z) 



+ 



4vr 2 
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29\ ixl 



5(1 -z) + 



8(1-2) 
772 + 27 An(l-z)' 



272 2 - 572 -2(772 + 27) In * 



1-2 



■3(i; 2 28; I- y - y ) : -) 



+ ^282 2 - 22,2 + 



107 18 
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8 4 82 



3s \ A 



m. 



m 2 z 2 + (1 - zf (z^ 17z _1 

11 /4 2 + v 4 + 8 4(1 + ^) 
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a cc (z,/j, F ) = 



327ra 2 /V ^ 2 



9s 



tc 1 2 

A m * 



, m 2 1 + (1 - zf 

ln^f + 2zln(l - 2) - 2 In 2 
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a C c(z,/j, F ) = 



167TO 2 (n 9 tc \ 2 

"9T- T m < 



, m? 1 + (1 - z) 2 n , 

ln-f i '- + 22 In 1 - 2) - 2ln2 

\i F z 



2 4 7 4 2 2 
--2 - 41n(l - 2) + 2hi2 1 h -ln(l - 2) In 2 + -z 2 

3 32 3 2 2 3 



Vcq(q)(z,V>F) = 



167ra 2 /«&.\ 2 

— mf 

9s A * 



, mf 1 + (1 - z) 2 n , 

ln^f i ^ + 22 In 1 - 2) - 2ln2 

/4 ^ 



4 2 2 

-4 ln(l - 2) + 2 In 2 + 3 + - ln(l - 2) - - In 2 - - 

2 22 



(19) 



Finally, we combine these finite a a b(z,fi F ) with the appropriate partonic distribution 
functions to arrive at the NLO cross sections for the direct top productions: 



a NLO = I dxidx 2 \ Gg/p 1 (x 1 ,i2 F )G c /p 2 (x 2 ,lJ,F) + (xi ^ x 2 ) a gc (z,n F ) 
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+G g/Pl (Xi, r l F )G g /p 2 (x 2 , fi F )Vgg(z, fI F ) + G c / Pl (xi, fi F )G c /p 2 (X 2 , [i F )Vcc(z, Vf) 

+ G c / Pl (xx, fi F )G 5 /p 2 (x 2 , Hf) + {xi <-> x 2 ) & c5 (z,fi F ) 

+ Y1 [ G e/Pi(Xl, fl F )G q (q ) /p 2 {x 2 , fJ, F ) + (SCI <-> 2C 2 ) ^cq(q){z,fi F ) >■ (20) 
9(5) J 

IV. NUMERICAL CALCULATION AND DISCUSSION 

In this section, we present the numerical results of the NLO QCD corrections to the direct 
top quark production via anomalous FCNC couplings at the Tevatron Run 2 (\/S = 2.0 TeV) 
and the LHC (y/S = 14 TeV). In our numerical calculations, we take the top quark mass 
m t = 178.0 GeV Q and the two-loopevolution of a a (fj, r ) Q witn ot,{M z ) = 0.118. 
Moreover, CTEQ6L (CTEQ6M) PDFs (ifij are used in the calculation of the LO (NLO) 
cross sections. As for the anomalous couplings, which appear in the expressions as quadratic 
factors, we choose -f- = 0.01 TeV . 

In Table HI we list the LO and NLO cross sections for the center value of the factorization 
and renormalization scale //p = fi r = m t . One can see that for gc — > t (gu — > t), the NLO 
predictions enhance the LO results 61% (59%) at the Tevatron Run 2, and the 40% (52%) 
at the LHC, respectively. 

We define R as the ratio of the LO (NLO) cross sections to their values at fi F = \i r = 
m t . Fig. El shows the ratio R as functions of factorization (renormalization) scales for sub- 
processes gc — > t and gu — > t at the Tevatron Run 2. For each subprocess, we show the 
scale dependence for three cases: (1) fi F and \i T varies simultaneously, (2) \i F = m t and 
only fi r varies, (3) fi r = m t and only fi F varies. We can find that, for the three cases, the 
NLO QCD corrections remarkably reduce the dependence of theoretical predictions on the 
factorization (renormalization) scales. For example, in the region m t /2 < /j, F = /i r < 2m t , 
the cross sections vary by ±18% at LO but by ±11% at NLO for gu — > t, and vary by ±15% 
at LO but by ±9% at NLO for gc -> t. 

Fig0 shows the ratio R as functions of factorization (renormalization) scales at the LHC. 
One can find that the scale dependence is reduced after considering the NLO QCD correc- 
tions in case (2) and (3). For case (1), although the NLO results show the weaker scale 
dependence when fj, F — [i r is large (> m t ), the NLO results do not reduce the scale depen- 
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subprocess 


LHC (LO) 


LHC (NLO) 


Tevatron Run 2 (LO) 


Tevatron Run 2 (NLO) 


gu -> t 


11069.8 


16817.8 


259.0 


412.6 


gc t 


1817.1 


2536.6 


17.6 


28.3 



TABLE I: The LO and NLO cross sections of direct top quark production via anomalous FCNC 
couplings at the LHC and Tevatron Run 2 (fb). Here -P- = 0.01 TeV -1 and \xf = = TfH- 



dence of theoretical predictions for fip = \x T < mt. The main origins to the scale dependence 
come from the terms proportional to 5(1 — z) in a gc (z, fip), and higher order effects are nec- 
essary to further reduce such scale dependence in the case of fip — n r at the LHC. Similar 



phenomena also can be found in other processes 
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V. CONCLUSIONS 



We have calculated the NLO QCD corrections to the cross sections for direct top quark 
productions induced by model-independent FCNC couplings at hadron colliders. The NLO 
results increase the experimental sensitivity to the anomalous couplings. Our results show 
that the NLO QCD corrections enhance the LO total cross sections at the Tevatron Run 
2 about 60% for both of k\ c and k\ u couplings, and enhance the LO total cross sections at 
the LHC about 40% for n 9 tc couplings and 50% for k\ u couplings, respectively. Moreover, 
the NLO QCD corrections vastly reduce the dependence of the total cross sections on the 
renormalization or factorization scale, which leads to increased confidence in predictions 
based on these results. 
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FIG. 2: Feynman diagrams of real gluon emission sub-processes ((ri) — (re)), gluon initial state 
sub-processes ((51) — (54)), of which (g±) and (52) have cross diagrams, and quark initial state 
sub-processes ((si) — (-S3)), of which (S2) has cross diagram. 
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FIG. 3: The ratio R as functions of factorization (renormalization) scales for sub-processes gc — > t 
and gu — > t at the Tevatron Run 2: (a) up quark initial state and (b) charm quark initial state. 
Here %^ = O.OlTeV" 1 . 
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FIG. 4: The ratio R as functions of factorization (renormalization) scales for sub-processes gc 
and gu — > t at the LHC: (a) up quark initial state and (b) charm quark initial state, f 
%4 = O.OlTeV- 1 . 
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